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Abstract 

As in many other locations in the world, honeybee colony losses and disorders have increased in Belgium. Some of the 
symptoms observed rest unspecific and their causes remain unknown. The present study aims to determine the role of both 
pesticide exposure and virus load on the appraisal of unexplained honeybee colony disorders in field conditions. From July 
201 1 to May 2012, 330 colonies were monitored. Honeybees, wax, beebread and honey samples were collected. Morbidity 
and mortality information provided by beekeepers, colony clinical visits and availability of analytical matrix were used to 
form 2 groups: healthy colonies and colonies with disorders (n = 29, n = 25, respectively). Disorders included: (1) dead 
colonies or colonies in which part of the colony appeared dead, or had disappeared; (2) weak colonies; (3) queen loss; (4) 
problems linked to brood and not related to any known disease. Five common viruses and 99 pesticides (41 fungicides, 39 
insecticides and synergist, 14 herbicides, 5 acaricides and metabolites) were quantified in the samples.The main symptoms 
observed in the group with disorders are linked to brood and queens. The viruses most frequently found are Black Queen 
Cell Virus, Sac Brood Virus, Deformed Wing Virus. No significant difference in virus load was observed between the two 
groups. Three acaricides, 5 insecticides and 13 fungicides were detected in the analysed samples. A significant correlation 
was found between the presence of fungicide residues and honeybee colony disorders. A significant positive link could also 
be established between the observation of disorder and the abundance of crop surface around the beehive. According to 
our results, the role of fungicides as a potential stressor for honeybee colonies should be further studied, either by their 
direct and/or indirect impacts on bees and bee colonies. 
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Introduction 

The evolution of pollinator populations has been the subject of 
an increasing number of studies, most of them showing worrying 
negative trends [1^]. Furthermore, beekeepers have long been 
notifying enhanced bee winter losses and disorders [5-8]. Bees 
contribute to ecosystem services and their decline thus threatens 
pollination of both wild and cultured plants, endangering 
biodiversity, food and fibre production [9,10]. This decline may 
also have an impact on the production of other goods with 
pharmacological uses (e.g. honey, propolis) [11,12], and scientific 
and technological inspiration (e.g. development of visual guided 
flight robotics) [13]. Bees are also part of our culture (e.g. culinary, 
hobby occupation, etc.), contributing to the dynamism of rural and 
urban areas [14] and providing a source of inspiration and well- 
being for many [15]. 

Belgium shares the trends observed worldwide both in terms of 
wild and reared pollinators [8,16], with enhanced winter mortality 
observed from 1999 [8,17,18]. However, apart from colony 
mortality, Walloon beekeepers have reported a number of 
imprecise symptoms: colony weakness, mainly in spring; fast 
renewal of young queens; rapid loss of individuals in the colony, 
mainly foragers, or slow loss of individuals in the colony. In many 



cases brood and food remains in the colony. Sometime a small 
cluster of bees with the queen survives [17]. Some beekeepers 
described unspecific brood abnormalities not characteristic to any 
known disease (Baudoin and Lequeux, pers. com.). Other studies 
describe similar mortality trends [8,19,20] as well as unspecific 
symptoms: increased colony mortality and/or weakening [21-23], 
queen failure [21,23-25] or honey yield reduction [21]. 

Multifarious factors have been proposed to provide a cause to 
such a phenomenon. Climate change is proposed as one of them, 
together with a decrease of genetic variability of the bee colonies, 
electromagnetic radiation, pathogens and parasites or the impact 
of intensive agricultural systems (nutritional lack, GM plants or 
pesticides) [26]. Previous studies developed in the region discard 
some of these factors as causes of bee losses, specifically Nosema 
spp and American Foulbrood [27]. However, pesticide residues 
and certain viruses were detected in bee colonies [28,27]. These 
elements were the most relevant to us. 

Coundess studies have shown lethal and sub-lethal effects of 
pesticides on bees [29-31]. Insecticides are often the most studied 
for obvious reasons. However other substances (fungicides, 
herbicides) deserve analysis for their specific toxicity, individually 
or in synergy with other substances [32,33] or pathogens [34-36], 
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or their extensive exposure given their large scale and/ or repeated 
use. Indeed, residue analyses show that these types of substance are 
found in the hives, even though the crops in which they were 
applied would suggest no bee exposure [37-40]. 

Beside pesticides, bee viruses are often mentioned as a putative 
cause of decline- of the colonies, or at least to reside among the 
presumed multi-factorial causes [41]. In Europe, at least 12 viruses 
infecting bees have been compiled [21]. Many honeybee viruses 
commonly occur in seemingly healthy populations that continue to 
run well. Viruses may remain latent and confined within certain 
cells or tissues with no active replication and no disruption of 
ceUular function. Likewise, they may be replicating at low level in 
permissive cells but in non-vital sites or in honeybee life stages that 
do not exhibit any symptoms or obvious pathology [21,42,43]. 
Nevertheless, two viruses, ABPV and DWV, are able of inducing 
serious disorders to honeybees and ha\'e been shown to cause -in 
association with V. destructor - winter losses in Germany [23]. 

In this paper, we first studied the relationship of in-hive viral 
prevalence as well as the presence of pesticide residues on 
honeybee colonies' health. As a next step, we focused on the 
relationship between the environment surrounding the monitored 
apiaries and the health condition of their colonies. 

Materials and Methods 

Field work - colony follow up 

A group of voluntary beekeepers were requested to participate 
in the study. All of them shared the following criteria: (1) have a 
minimum of 5 production colonies per apiary in June 20 1 1 ; (2) 
regularly follow up the- health status and development of their 
colonies; and (3) monitor the varroa infestation level and carry out 
officially recommended varroa treatments (treatment in July- 
August with veterinary medicaments based on th)^!^ and winter 
treatment with veterinary medicaments based on oxalic acid). A 
total of 330 colonies distributed among 66 apiaries (5 colonies/ 
apiary) in Walloon and Brussels regions (Belgium) were followed. 
In those apiaries composed of 5 or more colonies, the 5 well- 
developed colonies at the beginning of the study were selected. 

Colonies received three visits along the study, the first one from 
mid-July to mid- August 2011, the second one from mid- 
September to mid-October 2011 and the third one from March 
to April 2012. State official beekeeping technicians speciafised in 
bee health were trained in the framework of the project in order to 
minimise as far as possible the variability of the results due to 
handling and observations. For each of the visits, beekeepers were 
requested to fulfil a form asking information about the health 
history of the apiary and of the foUowed-up colonies, as well as 
about their beekeeping practices. Honeybee colony disorder 
symptoms were reported in the questionnaire. These include the 
following symptoms for which no causal agent could obviously be 
identified: 

(1) dead and disappeared colonies : (la) death of part or the 
whole colony, where dead bees can be found close to the hive 
or inside it. Beekeepers also describe the (lb) disappearance of 
part or the whole colony, leaving behind food reserves and 
brood, a phenomenon similar to the one described in North- 
American apiaries (Colony Collapse Disorder (CCD) [24]; 

(2) weak colonies : weakening of the colony, showing in occasions 
a slow development in spring under optimal conditions (e.g. 
optimal weather, low varroa pressure, etc) with as conse- 
quence the loss of the spring production, but in occasions 
showing abnormally small colonies or colonies with low 
activity; 



(3) queen loss : replacement of young queens sometimes leading 
to queen-less colonies or interruption of the egg-laying activity 
of the queen [17]; 

(4) problems linked to brood and not related to any known 
disease. 

The form also included requests about other symptoms typically 
linked to known diseases (e.g. diarrhoea, mummified larvae, 
varroa presence, presence of deformed wing bees, etc). Each of the 
questions requested further characterisation of the symptoms (e.g. 
population size, population dynamics, bees behaviour, etc.). 

In addition, a thorough clinical inspection was carried out for 
each colony in the mentioned period. Special attention was given 
to the strength of the colony in terms of bee numbers, brood 
surface and reserves content, the presence of the queen and the 
varroa infestation level. Finally, two different samples were 
collected before and after the winter: (1) in-hive bees (minimum 
of 100 bees); (2) a section of the frame containing beebread and 
honey (about 100 cm^). Samples were collected in hermetic plastic 
bags, cooled after collection and stored at — 20°C. 

Case choice - hierarchical sample clustering 

Information available from bee colonies comprised field 
observations, beekeeper answers to the questionnaire and results 
from the analyses of the samples of beekeeping matrices (honey, 
beebread, wax and bees). Colonies for which this information was 
available were considered for our analyses. Colonies with well 
identified problems (heavy varroa infestation, lack of food or 
drone-laying queens) were discarded. After that, two groups were 
made: a group with disorders and another with healthy colonies. 
In order to constitute these groups and to limit variations due to 
potential difierent bee management, a hierarchical classification 
was made. The criteria used were the amount of food stores before 
the winter, the year of colony creation, the subspecies and the age 
of the queen. This classification was made by using Ward 
aggregation method. This method allows to buUd group with the 
lower variation within the group and the higher variation between 
groups [44]. 

Virus analyses 

The viruses under investigation were the Black Queen Cell 
Virus (BQCV), the Chronic Bee Paralysis Virus (CBPV), the Acute 
Bee Paralysis Virus (ABP\'), the Deformed Wings Virus (DWV), 
and the Sac Brood Virus (SBV). Viral analyses were conducted on 
the honeybee workers collected before the winter, during the first 
and second periods July-August and September-October). The 
samples were analysed with a quantitative RT-PCR by the 
National Bee Unit laboratory. Food and Environment Research 
Agency (Sand Hutton, York, United-Kingdom). 

Total Nucleic acid (TNA) was extracted from 60 bees 
homogenised for 12 minutes with 20 ml GITC Lysis Buffer 
(5 M Guandine Thiocynate, 0.05 M Tris base, 0.02 M EDTA, 
pH 8.0) in a 30 ml bottie containing 3, 7/16" ball bearings. GITC 
Lysis buffer also contained 17.3 mM SDS buffer (173 mM 
Sodium dodecyl sulphate (SDS) in 100 ml MGW). The SDS 
buffer is added prior to use in warmed GITC Lysis Buffer. The 
homogenate was then incubated at 65°C for 40 minutes and then 
spun at 6189 g for 5 minutes. Polypropylene 96 -deep well plates 
(DT850301 Elkay Laboratory Products Ltd) were prepared as 
follows (1 well/ extract); plate A: 800 |J,1 extract, and 100 nl 
MagneSn™ beads (MD1441, Promega); plate B: 1 ml of GITC 
wash buffer (5 M Guandine Thiocynate, 0.05 M Tris base); plates 
C, D, E: 1 ml 70%v/v edianol (E/00665DF/17, Fisher Scientific); 
plate F: 1 ml 5 M Betaine (B2629, Sigma), plate G, 300 ul 1 xTE 
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Buffer (stock lOOX TE (EC-862 National Diagnostics). Plates were 
loaded onto the Kingfisher Flex and processed as follows: Plate A - 
Bind 5 mins (fast dual mix), Plate B - Wash 3 mins (fast dual mix), 
Plate C, D, E — Wash 2 mins (fast dual mix), collect beads at 1 min 
intervals, Plate F - Wash 20 sees (medium), without releasing 
beads, Plate G - Mix 1 min then incubate at 65°C for 5 mins with 
mixing. AH steps of the process are looped through twice. TNA 
was collected from plate G of each reaction and stored at — 80°C 
prior to use in real-time PGR. 

Reactions were set up in 96 well reaction plates using Absolute 
Blue QPCR ROX mix (AB-4139, Thermo Scientific) following the 
protocols supplied. 0.1 mM of 0.1 M Dithiothreitol (165680250, 
Arcos Organics) and 0.33 unit of MMLV (EP0441: Fermentas) 
were added to each reaction. The primers (Table 1) were all used 
at 400 nM and probes at 200 nM final concentration. Total 
nucleic acid (5 |J,1) was added to each reaction, giving a final 
reaction volume of 20 |xl. Plates were cycled for 48°C/30 min, 
95°C/10 min and 40 cycles of 60°C/1 min, 95°G/15 sec within 
the 7900 HT Sequence Detection System (Apphed Biosystems), 
using real time data collection. The results were recorded as the 
cycle threshold (Gt) or cycle number after which a significant 
accumulation of florescence over the baseline was observed; an 
average (of duplicate wells) Gt value below 40 was regarded as a 
positive result with a threshold ARn setting of 0.2. Given the 
absence of internal standard, we assumed that the extraction 
method had the same efficiency towards bee DNA/RNA and viral 
RNA. 

Virus primers and probes used for aU pathogens tested have 
been previously described in Chantawannakul et al., 2006 [45]. 
Additional APBV primers are described in Martin et al., 2012. 
[46]. 

Pesticide analyses 

Samples collected before the winter were sent to Eurofins 
Ghemiphar NV, Brugge, Belgium and analysed by SOFIA GmbH 
Ghemisches Labor fUr Softwareentwicklung und InteUigente 
Analytik, Berlin, Germany. Two multi-residues methods (SF146 
and SF150) were used searching for 99, 93 and 96 pesticides 
residues in wax (54 samples), beebread (108 samples) and honey 
(107 samples - one sample did not contain enough matrix) 
respectively (LOQ in Table 2). 

For SF146 method, 10 ml of water was added to 10 g of 
samples. Methanol was added and the preparation was mixed. 
The mixture was filtered and centrifuged. Next, for analysis by gas 
chromatography coupled with mass spectrometry (GC-MS), the 
filtrate was mixed with sodium chloride and ethyl acetate solution 



1:1 (v/v). All was dried with sodium sulfate and filtrated. This 
filtrate was then concentrated and analysed by GG-MS. For 
analysis with liquid chromatography coupled with tandem mass 
spectrometry (LG-MS/MS), 5 ml of sample were transferred in 
GhemElut column and eluted with dichloromethane. After 
concentration, 1 ml of water-methanol 1:1 (v/v) was added and 
this solution was analysed by LG-MS/MS. For the wax, some 
modifications were made. For GG-MS, sodium chloride was 
replaced by hexane to dissolve wax and the mixture was frozen 
until a precipitate appears. For LG-MS/MS, 0.2 g of sample were 
extracted with 30 ml of mixture of naphtha:acetonitrile 1:2 (v/v). 
1 0 ml of acetonitrile phase was isolated and concentrated before 
the analysis by LG-MS/MS. 

For SF150 method, 10 ml of water was added to 10 g of 
samples as well as methanol and hydrochloric acid. AU was mixed, 
filtrated and centrifuged. A solution of sodium chloride (20%) was 
added and 5 ml of sample were transferred in GhemElut column 
and eluted with dichloromethane. After concentration, 1 ml of 
water-methanol 1:1 (v/v) was added and this solution was analysed 
by LG-MS/MS. 

Descriptive analyses 

A descriptive analysis of virus prevalence and pesticide residues 
was first carried out. The relationship between disorder and these 
stressing factors was then assessed. For each virus type, a linear 
model (two way ANOVA) was used, with the cycle threshold value 
as dependent variable and the visit (first or second), group (with 
disorder or healthy) and their interaction as explanatory variables. 
For the pesticides residues, we used a similar model with the total 
number of residues as dependent variable and type (fungicide or 
insecticide/acaricide), group (with disorder or healthy) and their 
interaction as explanatory variable. Based on these models, we 
constructed a contrast matrix to test explicit post-hoc hypotheses 
(e.g. compare virus load between groups within each visit...). As 
several dependent variables had a strongly asymmetric distribution 
(i.e. non-normal), we computed all p-values by permutation 
(n = 1000). We applied a Bonferoni correction on the post-hoc tests 
p-values to take into account the multiplicity of the tests. 

Relationship between bee colony disorders and potential 
stressors or surrounding environment of the apiary 

We used two separate generalized linear mixed models with a 
binomial distribution and logit link function. Li both models, the 
"group" (with disorder or healthy) was used as dependent binary 
variable and the apiary was used as random effect (grouping 



Table 1. 


List of primers used for virus analyses. 






Target 


Primer Name 


Sequence (5 -3 ) 


BQCV 


BQCV 9195F 


GGT GCG GGA GAT GAT ATG GA 




BQCV 8265 R 


GCC GTC TGA GAT GCA TGA ATA C 




BQCV 821 7T 


TTT CCA TCT TTA TCG GTA CGC CGC C 


SBV 


SBV 31 IF 


AAG TTG GAG GCG CGy AAT TG 




SBV 380R 


CAA ATG TCT TCT TAC dAG AGG yAA GGA TTG 




SBV 331T (MGB) 


CGG AGT GGA AAG AT 


CPV 


CPV 304F 


TCT GGC TCT GTC TTC GCA AA 




CPV 371 R 


GAT ACC GTC GTC ACC CTC ATG 




CPV 325T 


TGC CCA CCA ATA GTT GGC AGT CTG C 


doi:l 0.1 371 /journal.pone.Ol 03073.t001 
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Table 2. Active ingredients or 
F = Fungicide; H = Herbicide; 1 = 


metabolite included in the multi-residue 
Insecticide; S = Synergist). 


analyses by beekeeping matrix (A 


= Acaricide; 




Active ingredient or metabolite* 


Class 


LOQ (mg/kg) 










Wax 


Beebread 


Honey 


2,4-D 


H 


0.1 


- 


0.01 


Abamectin 


1 


0.1 


0.1 


0.005 


Acetamiprid 


1 


0.1 


0.1 


0.005 


Aldicarb 


1 


0.1 


0.1 


0.005 


Alpha-cypermethrin 


1 


0.1 


0.1 


0.005 


Amitraz 


A 


0.1 


0.1 


0.01 


Azoxystrobin 


F 


0.1 


0.1 


0.005 


Bentazone 


H 


0.2 


- 


0.02 


Benthiavalicarb 


F 


0.1 


0.1 


0.005 


Beta-cyfluthrin 


1 


0.1 


0.1 


0.01 


Bifenthrin 


1 


0.1 


0.1 


0.005 


Boscalid 


F 


0.1 


0.1 


0.005 


Captan 


F 


0.1 


0.1 


0.01 


Carbaryl 


1 


0.1 


0.1 


0.005 


Chlorpyriphos 


1 


0.1 


0.1 


0.005 


Chlorpyriphos-methyl 


1 


0.1 


0.1 


0.005 


Chlorothalonil 


F 


0.1 


0.1 


0.005 


Clothianidin 


1 


0.1 


0.1 


0.01 


Coumaphos 


A 


0.1 


0.1 


0.05 


Coumaphos oxon * 


A 


0.1 


0.1 


- 


Coumaphos phenolic* 


A 


0.1 


0.1 


- 


Cyazofamid 


F 


0.2 


0.2 


0.02 


Cyfluthrin 


1 


0.1 


0.1 


0.01 


Cymoxanil 


F 


0.1 


0.1 


0.01 


Cypermethrin 


1 


0.1 


0.1 


0.005 


Cyproconazole 


F 


0.1 


0.1 


0.005 


Cyprodinil 


F 


0.1 


0.1 


0.01 


DDT 


1 


0.1 


0.1 


0.005 


Deltamethrine 


1 


0.1 


0.1 


0.005 


Dichlorprop-P 


H 


0.1 


- 


0.01 


Difenoconazole 


F 


0.1 


0.1 


0.005 


Diflubenzuron 


1 


0.2 


0.2 


0.02 


Dimethenamid-P 


H 


0.1 


0.1 


0.005 


Dimethoate 


1 


0.1 


0.1 


0.005 


Dimethotmorph 


F 


0.1 


0.1 


0.01 


Dimoxystrobin 


F 


0.1 


0.1 


0.005 


Epoxiconazole 


F 


0.1 


0.1 


0.005 


Esfenvalerate 


1 


0.1 


0.1 


0.005 


Ethofutmesate 


H 


0.5 


0.5 


0.05 


Famoxadone 


F 


0.1 


0.1 


0.01 


Fenhexamid 


F 


0.1 


0.1 


0.005 


Fenoxycarb 


F 


0.1 


0.1 


0.005 


Fenpropidin 


F 


0.1 


0.1 


0.01 


Fenpropimorph 


F 


0.1 


0.1 


0.005 


Fipronil 


1 


0.1 


0.1 


0.01 


Flonicamid 


1 


2 


2 


0.2 


Fluazinam 


F 


0.1 


0.1 


0.01 


Flufenacet 


H 


0.1 


0.1 


0.01 
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Table 2. Cont. 





Active ingredient or metabolite* 


Class 


LOQ (mg/ltg) 










Wax 


Beebread 


Honey 


Fluopicolide 


F 


0.1 


0.1 


0.01 


Flusilazole 


F 


0.1 


0.1 


0.005 


Heptenophos 


1 


0.1 


0.1 


0.01 


tmidacloprid 


1 


0.1 


0.1 


0.005 


Indoxacarb 


1 


0.1 


0.1 


0.005 


Iprodione 


F 


0.1 


0.1 


0.005 


Kresoxim-methy! 


1 


0.1 


0.1 


0.005 


Lambda-cyhalothrin 


1 


0.1 


0.1 


0.005 


Lindane 


1 


0.1 


0.1 


0.005 


Linuron 


H 


0.1 


0.1 


0.005 


MCPA 


H 


0.1 


_ 


0.01 


Mecoprop-P 


H 


0.1 


_ 


0.01 


Metalaxyl-M 


F 


0.1 


0.1 


0.005 


Metamitron 


H 


0.1 


0.1 


0.005 


Metconazole 


F 


0.2 


0.2 


0.02 


Metliiocarb 


1 


0.1 


0.1 


0.005 


Methoxyfenozide 


F 


0.1 


0.1 


0.005 


Metrlbuzin 


H 


0.1 


0.1 


0.005 


Penconazole 


F 


0.1 


0.1 


0.005 


Pendimethalin 


H 


0.1 


0.1 


0.005 


Permethrin 


1 


0.2 


0.2 


0.02 


Picoxystrobin 


F 


0.1 


0.1 


0.005 


Piperonyl butoxide 


S 


0.1 


0.1 


0.01 


Pirimicarb 


1 


0.1 


0.1 


0.005 


Prochloraz 


F 


0.1 


0.1 


0.005 


Propamocarb 


F 


0.1 


0.1 


0.005 


Propiconazole 


F 


0.1 


0.1 


0.01 


Propyzamide 


F 


0.1 


0.1 


0.01 


Pymetrozine 


1 


0.1 


0.1 


0.005 


Pyraclostrobin 


F 


0.1 


0.1 


0.005 


Pyrethrin 


1 


0.2 


0.2 


0.02 


Pyrimethanil 


F 


0.1 


0.1 


0.005 


Rotenone 


1 


0.1 


0.1 


0.01 


Spinosad 


1 


0.1 


0.1 




Spirodiclofen 


1 


0.1 


0.1 


0.01 


Spirotetramat 


F 


0.1 


0.1 


0.01 


Spiroxamine 


F 


0.1 


0.1 


0.005 


Sulcotrione 


H 


0.1 




0.01 


Tau-fluvalinate 


lA 


0.1 


0.1 


0.01 


Tebuconazole 


F 


0.1 


0.1 


0.005 


Tebufenozide 


1 


0.1 


0.1 


0.005 


Tefluthrine 


1 


0.1 


0.1 


0.005 


Terbuthylazine 


H 


0.1 


0.1 


0.01 


Tetraconazole 


F 


0.1 


0.1 


0.01 


Thiabendazole 


F 


0.1 


0.1 


0.005 


Thiacloprid 


1 


0.1 


0.1 


0.005 


Thiamethoxam 


1 


0.1 


0.1 


0.005 


Thiophanate-methyl 


F 


0.1 


0.1 


0.005 


Trifloxystrobin 


F 


0.1 


0.1 


0.01 
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Table 2. Cont. 



Active ingredient or metabolite* Class "-OQ (mg/kg) 

Wax Beebread Honey 

Zeta-cyperrrethrin I 0.1 0.1 0,005 

Zoxamide F 0.1 0.1 0.01 



doi:l 0.1 371/journal.pone.Ol 03073.1002 

factor) to take into account the non-independence between 
colonies from the same apiary and allow the intercept to vary 
between apiaries [47]. 

With the first model, we explored the relationship between the 
probability of disorder in a colony and potential stress factors, i.e. 
pesticides or viruses total load into the colonies. We used as 
explanatory variables the total number of fungicides residues, the 
total number of insecticides or acaricides (pooled together) 
residues, the total number of viruses detected for both visits and 
all first level interactions between these three explanatory 
variables. The insecticides and acaricides were pooled together 
because the most frequently found active ingredient, namely tau- 
fluvalinate, is authorized in Belgium for both purposes. 

With the second model, we investigated the relationship 
between the probability of disorder in a colony and the structure 
of the agricultural landscape around the beehive. The surfaces of 
all different kinds of agricultural soil occupancy in a circle with a 
radius of 1 500 m around each apiary were calculated according to 
the information provided by the farmers to the Walloon 
administration (Land Parcel Identification System). The different 
soil occupancy categories were pooled into two categories defined 
according to the potential frequency of plant protection product 
use: (1) grasslands (low pesticides use), (2) crops sensu lalo 
(potentially higher pesticides use), including major crops (mainly: 
cereals, potatoes, beet, oilseed rape, maize, flax, ...) but also 
surfaces dedicated to fruit or vegetables production, fodder 
production (mainly legumes) and to horticulture. The surfaces of 
these two land use options, grasslands and crops [sensu lata), were 
used as explanatory variables in the model. The 1500 m radius 
was chosen accordingly to the mean pollen and nectar collecting 
distance for the honeybee [48,49]. As the agricultural soil 
occupancy data was available only for WaUonia, we eliminated 
from the analysis the apiaries whose buffer had less than 70% of 
their surface within the administrative bovmdaries of WaUonia 
(n = 2). 

We used Likelihood Ratio (LR) Tests to evaluate the 
significance of the explanatory variables [47,50]. We respected 
the marginality rules, i.e. all main effects were tested after 
removing from the model the interactions in which they are 
involved [51]. 

AH analyses were performed with R (R Core Team 2013) and 
the mixed models were fitted with the package lme4 [52]. 

Results 

Description of honeybee colony disorders 

We gathered all sources of information, i.e. questionnaires, 
analytical matrices and information from clinical visits, for 173 
colonies. After the clustering, 54 colonies coming from 2 1 apiaries 
(Figure 1) were considered for the study: 25 presenting bee 
disorders and 29 not presenting them, defined as healthy group. 
Their data and samples were collected by 8 beekeeping 
technicians. 



In the group presenting disorders, 6 colonies were dead or 
unviable -with only a handful of bees remaining with the queen- at 
the spring visit (Table 3). Pre-wintering weakness, winter worker 
losses, late-season re-queening were reported in these cases. One 
colony was dead due to queen failure without re-queening in 
spring. Furthermore, a number of colonies were weak in spring. 
They were characterised by a low number of individuals and a 
slow development. In total, nine colonies lost their queens during 
the study, five of them leading to queenless colonies. These were 
considered as "disorder colonies" in the framework of the study. 
Finally, brood abnormalities not linked to known disease were 
observed in 1 0 colonies, most of them both before and after the 
winter. 

Symptoms of other diseases or bee parasite observation 
remained low within the selected colonies. Wax moths were 
observed in two cases, one in each group. Varroa destructor was 
present in all colonies of the study. All of them were treated 
according to official veterinary advice. Diarrhoea was observed in 
one of the colonies showing brood abnormalities. Anatomical 
modifications (i.e. deformed or undeveloped wings, small abdo- 
mens) were observed in 6 colonies, half of them being healthy and 
half of them showing disorders. No symptoms of any of the 
foulbroods were observed. 

Virus content 

One bees sample from the "disorder" group could not be 
analysed due to a lack of sufficient matrix. The three most 
abundantly found viruses were the Deformed Wings Virus (DWV), 
the Black Qiaeen Cell Virus (BQCV) and the Sac Brood Virus 
(SBV) (Figure 2). The Acute Bee Paralysis Virus (ABPV) and the 
Chronic Bee Paralysis Virus (CBPV) were found only in a very 




Figure 1. Spatial distribution of selected apiaries. 

doi:1 0.1 371/journal.pone.01 03073.g001 
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Table 3. Symptoms observed in the group with disorders. 





Symptoms of disorder 


Frequency 


Mortality 


2* 


Weakening 


3 


Queen failure 


5 


Brood problems 


9 


Mortality+Weakening 


2 


Mortality+Weakening+Queen failure 


2 


Mortality+Queen failure 1 


Queen failure+Brood problems 1 


Total 


25 



* One of these colonies not considered in the model due to a lack of viral results. 
doi:l 0.1 371 /journal.pone.Ol 03073.t003 
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Figure 2. Virus content according to the Cycle Threshold (CT) for the groups "with disorder" and "healthy". Boxplot of Cycle 
Threshold for the first and second visits (visit 1 - mid-July to mid-august and visit 2 - mid-September to mid-October) and the group with disorders 
(grey, n = 24 colonies) and the healthy one (white, n = 29 colonies). Deformed Wings Virus (DWV), Black Queen Cell Virus (BQCV), SacBrood Virus (SBV), 
Acute Bee Paralysis Virus (ABPV), Chronic Bee Paralysis Virus (CBPV). In red, mean with confidence interval estimated by boostrap method. Note: CT 
values below 40 were regarded as positive results and the lowest CT values correspond to the higher virus contents. 
doi:1 0.1 371/journal.pone.01 03073.g002 
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1 2 

Visit 

Figure 3. Average number of different viruses per colony. Data 
sliown for the first and second visits (visit 1 - mid-July to mid-august 
and visit 2 - mid-September to mid-October) and the group with 
disorders (grey, n = 24 colonies) and the healthy one (white, n = 29 
colonies). Whiskers show the standard error (SE). 
doi:1 0.1 371 /journal.pone.01 03073.g003 

limited number of samples and did not allowed particular 
statistical analyses on these two viruses. 

There was no significant difierence in viral content between the 
group with disorders or the healthy one, and independently of the 
visit for any of the three most abundant viruses (2-way ANOVA 
tested by permutations, DWV: p = 0.731, BQCV: p = 0.373, SBV: 
p = 0.54). We observed a decrease of virus content from summer to 
fall 2011, independently of the group, for BQCV (p<0.001) and 
SBV (p<0.001) but not for DWV (p = 0.271). For BQCV only, we 
observed a significant (p = 0.036) group xvisit interaction indicat- 
ing that the decrease of the virus abundance from visit 1 to visit 2 
was higher in the group "with disorder" than in the group 
"healthy". Overall, the number of virus changed between visit 1 
and visit 2 and a significant decrease was observed in the group 
with disorders (p<0.001) (Figure 3). 

Pesticide analysis 

172 agrochemical residues of 23 different active ingredients 
were detected in 94 out of 269 samples. Wax was the most 
contaminated matrbc: 109 residues of 1.5 different active ingredi- 
ents; while 39 and 24 residues of 10 and 8 substances were 
detected in beebread and honey, respectively (Figures 4 & 5). 




b. 



Figure 4. Proportion of samples containing residues of 
acaricides/insecticides in the different beeiceeping matrices 
(honey, beebread and wax). Data shown for the group with 
disorders (grey, n = 25 colonies) and the healthy one (white, n = 29 
colonies). 

doi:10.1 371/journal.pone.01 03073.g004 

Residue levels contained in wax and beebread were higher (0.21— 
3.1 mg/kg) than those in honey (0.001-0.058 mg/kg). 

For the subsequent statistical analysis, the results obtained for 
insecticide and acaricides residues were pooled together because 
the most frequently found active ingredient, tau-fluvalinate 
(n = 46), can be used in Belgium as an acaricide against varroa 
mite and as an insecticide to control Meligethes aeneus in rape. 
The second most frequently found residue was the coumaphos 
(n = 35), followed by two fungicides, boscalid (n = 33), iprodione 
(n=13) (Table 4). Some residues of neonicotinoid insecticides, 
synergist and herbicide were also detected: thiacloprid (n = 3), 
piperonyl butoxide (n = 6), terbuthylazine (n=l). The highest 
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Figure 5. Proportion of samples containing residues of fungicides in the different beel<eeping matrices (Kioney, beebread and wax). 

Data shown for the group with disorders (grey, n = 25 colonies) and the healthy one (white, n = 29 colonies). 
doi:10.1371/journal.pone.0103073.g005 



residue level concerned captan with 3.1 and 1.9 mg/kg in the wax 
and bee bread of the same colony, respectively. Despite of being 
the most frequendy found in matrices, tau-fluvalinate and 
coumaphos residues never exceeded 0.71 and 0.58 mg/kg, 
respectively. Boscalid, the most commonly found fungicide, ranged 
from 0.005 to 1.3 mg/kg. 

There was a significant difference (p = 0.01) in terms of number 
of fungicide substances found between the "disorder" group 
(mean = 2.0) and the "healthy" one (mean = 0.7) as presented in 
Figure 6. The total number of insecticides/acaricides residues was 
slightly higher in the disorder group (mean = 2.1) than in the 
healthy group (mean = 1 .64) but this difference is not statistically 
significant (p = 0.79). The most frequent fungicides in the group 
with disorders were boscalid and iprodione, detected in the three 
investigated matrices (Figure 5). However, for single active 



substances, no significant difference of number of residues was 
observed between the groups "with disorder" and "healthy". 

Probability of disorders in relation to potential stressors 

There is clearly a significant positive relationship between the 
probability of a colony showing disorders and the total number of 
fungicides (Figwe 7, LR= 7.128, df = 1, p = 0.008). The estimated 
probability for a colony to be in the "with disorder" group is 0.26 
without fungicides residues, 0.60 with 2 fungicides residues and 
0.88 with 4 fungicides residues when the insecticides-acaricides 
residues number and virus load are fixed to their observed mean 
value (Figure 7). None of the other variables (total number of 
viruses, total number of insecticides-acaricides) and none of the 
first level interactions seem to have any explanatory power on the 
probability of disorders in a colony (Table 5). No direct link could 
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1 1 

Fungicides InsecticidesAcaricides 

Figure 6. Average number of residues per colony. Data shown 
for the group with disorders (grey, n = 25 colonies) and the healthy ones 
(white, n = 29 colonies). Whiskers show the standard error (SE). 
doi:10.1371/journal.pone.0103073.g006 

be established between bee colony disorders and the amount (in 
|ig/kg) of residues found in the matrices. 

Probability of disorders and crop/grassland surface 

Our data clearly show a significant increase in disorder 
probability with the increase of crop surfaces {sensu lato, i.e. 
including fruit, vegetables, fodder production and horticulture) in 
the surrounding of the apiary (Figiare 8, LR = 8.052, df= 1, 
p = 0.0045). The predicted probability of disorders is close to 0.1 
for a crop surface of 0 ha in the radius of 1500 m and increases up 
to 0.8 for a surface of 500 ha of crops when fixing the grasslands 
surface to its observed mean. 

On the contrary, the probability of disorders strongly decreases 
when the grassland surfaces increases (Figure 8, LR= 14.527, 
df= 1, p = 0.0001) after controlling for the crop surface. The 
predicted probability of disorders is close to 1 for a grassland 
surface of 0 ha in a radius of 1500 m and drops to ~0.1 for a 
surface of 150 ha of grasslands when fixing the crops surface to its 
observed mean. 

Very similar results were found with a 3000 m buffer around 
the apiary (results not shown). 




012345678 
Number of fungicide residues 

Figure 7. Probability of honeybee colony disorders depending 
on the number of fungicide residues detected. Model based on 
averaged coefficients and median value both for the number of 
insecticides-acaricides residues and total number of virus detected for 
both visits (n = 53 colonies). 
doi:1 0.1 371/journal.pone.01 03073.g007 

Discussion 

The bee disorders observed after the winter in the Walloon 
region happened despite of normal climatic conditions. The 
autumn of 201 1 was dry (140.4 mm rainfall water from October to 
December, average being 2 19.9 mm), sunny and warm for Belgian 
conditions (12.4°C on average, which is normally 10.9°C for this 
period), followed by average winter and spring 2012 in terms of 
temperature, rainfall, with the exception of a cold wave lasting 
twelve days in February [53]. These conditions would not induce, 
a priori, a risk for honeybee colonies. Furthermore, based on 
studies carried out on the palynological diversity of the pollens 
collected by these colonies before the winter (Table SI), the 
hypothesis of nutritional lack is unlikely. AH samples analysed 
contained at least 8 different botanic sources of pollen, including 
rich protein content as oilseed rape, ivy and Phacelia. 

Viruses infections has often been mentioned as a source of stress 
for honeybee colonies [7]. However no significant difference in 
quantity or occurrence was observed between the healthy and the 
"with disorder" groups. DWV is one of the most commonly 
detected virus mA. mellifera. The prevalence of this virus is even 
more important at colonies infested by Varroa destructor, a well- 
established vector of this virus [7] . In accordance with the present 
study, further studies run on Belgian apiaries show DWV, BQCV 
and SBV as the most frequently found viruses [54] . Unlike DWV, 
we observed that the amount of SBV and BQCV has dropped 
significantly between the first and second visits for each group. In 
the case of BQCV, this is of no surprise as the cycle of incidence of 
this virus has been shown to increase in late winter, with a peak in 
May or June followed by a rapid decline [55]. The observed 
decrease of SBV has also been reported in other studies [42] in 
which the authors suggested that bees could develop a molecular 
defensive mechanism to reduce virus multiplication, or that the 
change in bee susceptibility to the virus could result from 
environmental conditions such as the quality of pollen. 

We cannot prove any causal relationship between any of these 
viruses analysed and bee disorders, nor does the interaction 
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Table 5. Analysis of deviance table for generalized linear binomial models describing the relationships between the colony 
disorder probability and three variables: the total number of (1) fungicide residues (totfungicides), (2) insecticide-acaricide residues 
(totinsaca) and (3) virus detected for both visits (totvirus). 



LR df pOChisq) 



totfungicides 


7.13 


1 0.008 


totinsaca 


0.005 


1 0.943 


totvirus 


0.136 


1 0.712 


totfungicides:totvirus 


2.222 


1 0.136 


totinsaca:totvirus 


0.975 


1 0.323 


totfungicides:totinsaca 


0.901 


1 0.342 



LR = likelihood ratio. 

doi:l 0.1 371 /journal.pone.Ol 03073.t005 



between these two factors. Cox-Foster et al., 2007 [41] found no 
clear correlation between a variety of pathogens, including 
Nosenia spp., DWV, CBPV, ABPV, BQCV, Mellisococcm pluton 
and Paenihacillus larvae ssp and CCD. No specific spore counts 
were carried out in our study. However, no signs of nosemosis or 
foulbrood could be linked to colonies presenting bee disorders. It is 
noteworthy that a monitoring run at Belgian level found Nosema 
spp spores in one out of four colonies [56] . Nevertheless, Cox- 
Foster et al., 2007 [41], show a positive correlation between lAPV 
and CCD, which a priori would not be relevant in our conditions 
given the low prevalence of lAPV in this country [27,54]. A recent 
publication [36] shows a positive correlation between the presence 
of fungicides in pollen loads and the presence of spores of Nosema 
cei'anae. We do not exclude the potential involvement of Nosema 
spp. in the case of bee disorders. However, in the framework of our 
study, A^. ceranae seems either to play a role in the development of 
this weakening, while remaining asymptomatic, or not to play a 
decisive role in it. 

When considering pesticide residues, no direct link could be 
estabhshed between bee colony disorders and the amount (in |a.g/ 
kg) of residues found in the matrices. Neither could we identify any 
specific molecule as cause of bee disorders. Nevertheless, the study 
of the residue load of pesticides, specifically fungicides, opens new 
avenues for a better understanding of honeybee colony disorders. 

Even if insecticides/acaricides were the most numerous residues 
detected in hives mainly in wax, no significant difference in the 



number of accumulated residues was observed between colonies 
with disorders and the healthy ones. Indeed the two most 
abundant active ingredients, tau-fluvalinate and coumaphos, came 
most probably from varroa control measures even if tau-fluvalinate 
is used as an insecticide (Mavrik 2F) against Meligethes aeneus in 
oilseed rape. These active ingredients seem to be a frequent 
outcome of residue analyses studies [57,58,39,38,40,59]. Syner- 
gistic effects have been shown between acaricides and other 
molecules [33,60]. Nevertheless, our modelling did not suggest any 
synergistic effects in the appearance of bee disorders occurring 
between residues of fungicides and insecticides-acaricides. Resi- 
dues of synergist, piperonyl butoxide, were also detected that 
indicates a prior exposure to synthetic or natural pyrethroids even 
though residues were not found in the analysed matrix, probably 
due to a fast degradation of these active ingredients [61]. A non- 
authorised active ingredient was also detected: carbaryl, forbidden 
in Belgium since 2007, indicating an illegal agricultural or 
gardening use. 

The only neonicotinoid found, thiacloprid, was detected in 
honey during the sampling period of July-August in two colonies. 
However, the limits of detection achieved in our study do not allow 
to determine the presence of neonicotinoids -with the exception of 
acetamiprid- or fipronil at levels in the range of the acute toxicity 
of these substances (30-40 |J,g/Kg). Residues of these substances 
could be present at lower levels and thus an exposure to these 
substances cannot be excluded. 
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Figure 8. Probability of honeybee colony disorders depending on the apiary's environment. Consideration of crop surface vs. grassland 
surface in a radius of 1500 m around the apiary. Crops include fruit, vegetables, fodder production and horticulture. For each graph, the value of the 
variable not displayed is fixed to its observed mean (n = 53 colonies). 
doi:10.1371/journal.pone.0103073.g008 
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According to our results, the number of fungicide residues seems 
to plays a role in the appearance of honeybee colony disorder. 

Significandy, more fungicides residues were detected in colonies 
with disorders than in the healthy ones. Mainly four active 
ingredients, frequently used in plant protection, were detected: 
boscaUd, iprodione, pyrimethanil and cyprodinil. Fungicides are 
often considered safe for honeybees based on their acute toxicity. 
However some direct toxic effects on bees either by the mother 
compounds or their metabolites have been reported in the past 
[62]. For example, boscalid shows low acute toxicity to bees [63], 
despite the fact that beekeepers in the USA have already reported 
losses and adverse effects on bee brood development related to the 
foliar application of this systemic active ingredient [64]. Incidents 
were often linked to the use of a co-formulation boscalid - 
pyraclostrobin. However, both molecules appeared above detec- 
tion levels in only two colonies of our study and could not explain 
the general trend. Interestingly, one metabolite of boscalid, the 2- 
chloronicotinic acid, is similar to the 6-chloroiii(:otiiii(: acid, a 
metabolite of imidacloprid. The latter has proved to be lethal to 
bees at low concentrations (0. 1 Hg/L) following chronic exposure 
[6,5] . Further research should be carried out in order to clarify the 
mechanism of bee toxicity of boscalid. Indeed, boscalid has proved 
toxic to other aquatic invertebrates, r(;ducing daphnid fecundity 
and Chironomid emergence [64]. In addition, synergism with 
other active ingredients like insecticides are possible and increase 
the toxicity for honeybee [66,32,67]. Two other fungicides were 
also detected at very high levels in wax and beebread: iprodione 
and captan. The former is known for its synergistic effects in 
collaboration with insecticidal compounds [33]. The latter has 
been shown to induce effects on growth and development of larval 
honeybees [68,69]. Chronic and larval toxicity studies would be 
interesting at this stage in order to evaluate possible direct toxic 
effects on bee individuals and their behaviour. Indeed, a recent 
study showed increased larval mortality following chronic expo- 
sure to tau-fluvalinate, coumaphos, chlorpyriphos and chlorotalo- 
nyl or some of their combinations [70]. AU these substances were 
found in our study. Increased mortality rates in the fall may 
compromise the size and age structure of the wintering cluster, 
which could lead to winter losses. 

The indirect effects of fungicides on bees or on bee colonies are 
relatively little-known to date. 

Fungicides may have an impact on the colony by modifying the 
existing microflora present in the food stores or in the bee 
intestinal tract [71]. Studies have already shown the possible 
modification of microbial composition both at beebread level 
[62,72] and at intestinal level [73]. This modification in the 
composition of microbiota may lead to dysbiosis [74]. The impact 
that such an unbalance in the bee gut microflora may have on bee 
health has already been considered. The link between the 
unspecific symptoms observed in our study and a possible 
microbial alteration could be subject of further research. 

In parallel, the potential impact of microbial modification on 
digestibility and availability of nutrients should be a target for 
further research. Indeed, the content of essential amino acids 
might be altered when beebread is contaminated with fungicides 
(DeGrandi-Hoflfman, 2013, pers com.). Given the importance of 
nutrition, especially pollen, in the good development of the colony 
[75] alterations in composition or lack of essential nutrients would 
put the homeostasis of the colony at stake. Some studies have 
already shown the impact of nutritional lack on bee development 
and health [76]. Provided that poUen is the unique source of 
amino acids for honeybees, royal jelly production could also be 
affected [77,78] with unexpected potential consequences for its 
main consumers, larvae and the queen. A poor nutrition of the 



queen, could have as a result an impact in its activity. Likewise, a 
poor nutrition of the larvae has been shown to impact their 
development [76]. As a result, the presence of fungicides on 
beebread and honey may have both a direct effect on their health, 
but also an indirect one on the colony development. 

Fungicides are widely used in agriculture and are broadly 
present in bee matrices, sometimes at high concentrations at levels 
of mg/kg [79,80,38,59]. Boscalid, cyprodinil, iprodione are used 
to control a broad range of fungal pathogens including Botrytis 
spp., AUernaria spp. and Sclerotinia spp. on a wide range of crops 
including fruits, vegetables and ornamentals. Pyrimethanil is more 
specifically used to control grey mould on fruits, vegetables and 
ornamentals, and leaf scab on apple trees [81]. These active 
ingredients were already reported in bee matrices in Europe and in 
the USA [23,39,62]. Their frequent presence in bee matrices 
might be an indication of chemical plant protection intc-nsity in all 
agricultural landscapes. Fungicides could also be markers of 
exposure to other active ingredients with higher toxicity to bees. 
Mixes of products like fungicides — insecticides are often applied to 
reduce the number of spray applications. As a result, other 
pesticides often used in combination ^\ ith the fungicides found or 
applied in the same crops could have been at the origin of the 
effects observed in this study. However, the sensibility of the 
residue analyses used in our study might explain the lack of 
detection of such components. Further intensive monitoring and a 
thorough record of the agricultural practices on pesticide 
application in tank mixes would help clarifying this alternative 
explanation. 

Factors different from fiingicides are most likely involved in the 
development of bee disorders. According to our model, the 
disorder probability was not absent when fungicides residues were 
not detected in presence of insecticide residues and virus. Non 
identified pathogens, chemicals or factors of different nature could 
be operating as silent stressors. Ravoet and colleagues (2013) [82] 
reported the presence of new pathogens in Belgium that were not 
taken into consideration in our analyses (i.e. Crilhidia mellificae or 
the Lake Sinai Virus (LSV)). Other stressing factors could be also 
linked to the intensive agriculture. In fact, we observed an increase 
of colony disorders in the area with high density of crops in 
comparison with areas with grassland. Some studies and beekeeper 
claims go in line with our outcome regarding the concentration of 
bee problems in areas with intensive agriculture [83-86], 
Furthermore, the same negative trends on pollinators and 
biodiversity in areas with intensive agriculture have already been 
described as the result of as increased pesticide use, decreased 
landscape heterogeneity, loss and fragmentation of natural habitat 
[87,88]. 

In conclusion, the five virus studied (ABPV, CBPV, QBCV, 
SBV, DWV) do not seem determinant in the appearance of bee 
disorders in our study. These disorders seem clearly linked to the 
environment of the apiaries and were observed mainly in 
agricultural crop areas. We observed also that the number of 
fungicide residue's appears as the main potential stress factor linked 
to bee disorder. However other stressing factors could be acting or 
interacting at the same time: insecticides exposure, a lack of amino 
acids and oligo-elements, etc. Our results open new avenues for 
future research in order to better understand the side effects of 
fungicides on the bee colony and rjuestions the sustainability of 
intensive agriculture model and its impact on bees. Specific 
toxicological studies on both adult bees and larvae would be 
recommendable in order to better characterise the toxicity of 
fungicides. 



PLOS ONE I www.plosone.org 



14 



July 2014 I Volume 9 | Issue 7 | e103073 



Honeybee Colony Disorder in Crop Areas 



Supporting Information 

Table SI Diversity of pollens collected in the apiaries 
before the winter. 

(DOC) 

Acknowledgments 

We would like to salute and to thank the beekeepers and beekeeping 
technicians for their dedication and active contribution to our study. 
Likewise, we thank Gloria DiGrandi-Hoffman, Gerard Arnold, Szaniszlo 
Szokc and Marlin Derminc for their crilical and construclix'C rc\ic\v of iho 
article. We would like also to thank the Direction Gencralc Opcrationnelle 

References 

1. Biesmeijer JC, Roberts SPM, Reemer M, OhlemuUcr R, Edwards M, et al. 
(2006) Parallel declines in pollinators and insect-poUinated plants in Britain and 
the Netherlands. Science 313: 351-354. 

2. Potts SG, Biesmeijer JC, Kremen C, Neumann P, Schweiger O, et al. (2010) 
Global pollinator declines: trends, impacts and drivers. Trends Ecol Evol 25: 
345-353. 

3. Goulson D, Lvc QC, Darvill B (2008) Decline and conservation of bumble bees. 
Annu Rev Entomol 53: 191-208. 

4. Car\'alhriro LC, Kiinin WE. Kcil P. Ag-uirrc-CiiiliciTr/ j. Ellis \VN, ct al. (2013) 
Species richness declines and liiotic homogenisation ha\'e slowed down for NVV- 
European pollinators and plants. Ecol Ertt 16: 870-878. doi:10.1 1 1 1/ele. 12121 . 

5. Kluser S, Peduzzi P (2007) (iiobal pollinator decline: a literature review. UNEP/ 
GRID- Europe. © UNEP 2007. Switzerland. 

6. Maxim L, van der Sluijs JP (2010) Expert explanations of honeybee losses in 
areas of extensive agriculture in France: Gaucho ® compared with other 
supposed causal factors. Environ Res Lett 5: 014006. doi:10.1088/1748-9326/ 
5/1/014006. 

7. vanEngelsdorp D, Meixner MD (2010) A historical review of managed honev 
bee populations in Europe and the United States and the factors that mav afieet 
them. J Invertebr Pathol 103: S80-S95. 

8. Van der Zee R, Pisa E, Andonov S, Brodsehneider R, Charriere J-D, et al. 
(2012; Managed honey bee colony losses in Canada, China, Europe, Israel and 
Turkey, for the wmtcrs of 2008-9 and 1009-10. J Apie Res Bee World 51: 100- 
114. 

9. Klein A-M, Miiller C, Hoehn P, Kremen C (2009) Understanding the role of 
species richness for crop pollination services. Biodiversity, ecosystem function 
and human weUbeing. New York. 195-208. 

10. Lautenbaeh S, Seppelt R, Liebseher J, Dormann OF (2012) Spatial and 
temporal trends of global pollination benefit. PloS One 7: e35954. 

1 1. Banskota AH, 4'ezuka Y, Kadota S (2001) Recent progress in pharmacological 
research of propolis. Phytother Res 15: 561- 571. doi: 10. 1002/ptr. 1029. 

12. JuU AB, Rodgers A, Walker N (2008) Honey as a topical treatment for wounds. 
Cochrane Database of Systematic Reviews. Chichester, UK: John WUey & Sons, 
Ltd. 1-43. 

13. Srinixasan MV (2011) Honeybees as a Model for the Study of Visually Guided 
Flight, Navigation, and Biologically Inspired Robotics. Physiol Rev 91: 413—460. 
doi:10.1152/physrev.00005.2010. 

14. Bradber N (1990) Beekeeping in rural development. International Bee Research 
Association. 15 pp. 

15. UNEP (2005) Ecosystems and Human WcU-Bcing: Our Human Planet: 
Summary for Decision Makers. Island Press. Washington: Island Press. 

16. Rasmont P, Pauly A, Terzo M, Patiny S, Miehez D, et al. (2005) The survey of 
w'ild bees (Hymenoptera, Apoidea) in Belgium and Prance. Eood Agrie Organ 
Rome: 18. 

17. Lefevbre M, Bruneau E (2004) Etat des lieux du phcnomene de deperissement 
des ruchers en Region wallonne. Louvain la Neuve. 

18. Nguyen BK, MignonJ, Laget D, de GraafD, Jacobs F, et al. (2010) Honey bee 
colony losses in Belgium during the 2008-9 winter. J Apic Res 49: 337—339. 
doi: 10.3896/IBRA. 1 .49.4.07. 

19. Neumann P, Carreek NE (2010) Honev bee colony losses. J Apic Res 49: 1—6. 

20. Vanengelsdorp D, Hayes J, L'nderwood RM, Caron D, Pettis J (201 1) A survey 
of managed honey bee eolonv losses in the USA, fall 2009 to winter 2010. J Apic 
Res 50: 1-10. 

21. Aubcrt M, Faucon J-P, Chauzat M-P (2008) Eneiuete prospective multi- 
factorielle: influence des agents microbiens et parasitaires, et des residus de 
pesticides sur le devenir de colonies d'abeilles domesUqucs cn conditions 
naturelles. AFSSA. 

22. Higes M, Martin-Hernandez R, Garrido-Bail6n E, Gonzalez-Porto AV, Garcia- 
Palencia P, et al. (2009) Honeybee colony collapse due to Nosema ceranae in 
professional apiaries. Environ Microbiol Rep 1: 110-113. doi:10.1 1 1 l/j.l758- 
2229.2009.00014.x. 

23. Genersch E, von der Ohe W, Kaatz H, Schroeder A, Otten C, et al. (2010) The 
German bee monitoring project: a long term study to understand periodically 
high winter losses of honey bee colonies. Apidologie 41: 332-352. 



Agriculture, Rcssourccs Naiurcllcs ct Environnement (DG03), Dcpartc- 
ment des Aides (D4), Direction des Surfaces (D42), Service 42/3 - LPIS 
(Land Parcel Identification System) - Service Public de Wallonie for the 
spatial data. We would also like to thank the Molecular Technology' Unit at 
The Food and Environment Research Agency, Sand Hutton, York, Y041 
ILZ for Bee Pathogen Screening Services. Finally, we would like to thank 
the reviewers for their helpful comments. 

Author Contributions 

Conceived and designed the experiments: NS LH EB. Performed the 
experiments: NS LH CM. Analyzed the data: GSM LH. Wrote the paper: 
LH NS LAM GSM. 



24. vanEngelsdorp D, Speybrocck N, Evans JD, Nguyen BK, Mullin C, et al. (2010) 
Weighing Risk Factors Associated With Bee Colony Collapse Disorder by 
Classification and Regression Tree Analysis. J Econ Entomol 103: 1517-1523. 
doi:10.1603/EC09429. 

25. Pettis JS (2013) The Role Of Pesticides In Queen Health & Sperm Viability. 
Apimondia 2014, Kiev. 

26. Farooqui T (2013) A potential link among biogenic amines-based pesticides, 
learning and memory, and colony collapse disorder: A unique hypothesis. 
Neuroehem Int 62: 122-136. doi:10.1016/j.neuint.2012.09.020. 

27. Nguyen BK, Ribiere M, Vanengelsdorp D, Snoeck C, Sacgerman C, ct al. 
(20 1 1) Effects of honey bee virus prevalence, Varroa destructor load and queen 
condition on honey bee colony survival over the winter in Belgium. J Apic Res 
50: 195-202. 

28. Nguyen BK, Saegerman C, Pirard C, MignonJ, Widart J, et al. (2009) Does 
imidacloprid seed-treated maize have an impact on honey bee mortality? J Econ 
Entomol 102: 616-623. 

29. Desneux N, Deeourtye A, Delpueeh J-M (2007) The Sublethal Effects of 
Pesticides on Benelleial Arthropods. Annu Rev Entomol 52: 81—106. 
doi: 10.11 46/annurev.cnto.52. 1 1 0405.09 1 440. 

30. Belzunces LP, Tehamitchian S, BrunetJ-L (2012) Neural effects of insecticides in 
tiie honey bee. Apidologie 43: 348-370. doi:10.1007/sl3592-012-0134-0. 

31. Van der Sluijs JP, Simon-Delso N, Goulson D, Maxim L, Bonmatin J-M, et al. 
(2013) Neonicotinoids, bee disorders and the sustainability of pollinator services. 
Curr Opin Environ Sustain 5: 293-305. doi:10.1016/j.eosust.2013.05.007. 

32. Pilling ED, Jepson PC (1993) Synergism between EBI fungicides and a 
pyrethroid insecticide in the honevbee (Apis mellifera). Pestic Sci 39: 293—297. 

33. Fischer R, Waehendorlf-Xeumann U (2003) Active ingredient combinations 
with insecticidal, fungicidal and acaricidal properties. 

34. Alaux C, Brunet J-L, Dussaubat C, Mondet F, Tchamitchan S, et al. (2010) 
Interactions between Nosema microspores and a neonicotinoid weaken 
honeybees {Apis mellifera ). Environ Microbiol 12: 774—782. doi:10.1111/ 
j.l462-2920.2009.02123.x. 

35. Vidau C, Gonzalez-Polo RA, Niso-Santano M, Gomez-Sanchez R, Bravo-San 
Pedro JM, et al. (2011) Fipronil is a powerful uncoupler of oxidative 
phosphorylation that triggers apoptosis in human neuronal cell line SHSY5Y. 
NeuroToxieology 32: 935-943. doi: 10. 1016/j.neuro.201 1.04.006. 

36. Pettis JS, Lichtcnberg EM, Andree M, Stitzinger J, Rose R (2013) Crop 
pollination exposes honey bees to pesticides which alters their susceptibility to 
the gut pathogen Nosema ceranae. PloS One 8: e70182. 

37. Bogdanov S (2005) ContamiuEints of bee products. Apidologie 37: 1—18. 
doi: 1 0. 1 05 1 /apido:2005043 . 

38. Johnson RM, Ellis MD, Mullin CA, Frazier M (2010) Pesticides and honey bee 
toxicity - USA. Apidologie 41: 312-331. doi:10.1051/apido/2010018. 

39. Mullin CA, Frazier M, Frazier JL, Ashcraft S, Simonds R, et al. (2010) High 
Levels of Miticides and Agrochemieals in Xorth American Apiaries: Implications 
for Honey Bee Health^ PLoS ONE 5: e9754. doi:10.1371/journal.pone. 
0009754. 

40. Chauzat M-P, Martel A-C, Cougoule N, Porta P, Lachaize J, et al. (2011) An 
assessment of honeybee colony matrices, Apis mellifera (Hymenoptera: Apidae) 
to monitor pesticide presence in continental France. Environ Toxicol Chem 30: 
103-111. doi:10.1002/etc.361. 

41. Cox-Foster DL, Conlan S, Holmes EC, Palacios G, Evans JD, et al. (2007) A 
metagenomic survey of microbes in honev bee colony collapse disorder. Science 
318: 283-287. 

42. Tentcheva D, Gauthier L, Zappulla N, Dainat B, Cousserans F, et al. (2004) 
Prevalence and seasonal variations of six bee viruses in Apis mellifera L. and 
Varroa destructor mite populations in France. Appl Environ Microbiol 70: 
7185-7191. 

43. Gauthier L, Tentcheva D, Tournaire M, Dainat B, Cousserans F, et al. (2007) 
Viral load estimation in asymptomatic honey bee colonies using the quantitative 
RT-PCR technique. Apidologie 38: 426-^35. 

44. Zuur AF, leno EN, Smith GM (2007) The Analysis of Ecological Data. Springcr- 
Verlag. 

45. Chantawannakul P, Ward L, Boonham N, Brown M (2006) A Scientific Note on 
the Detection of Honeybee Viruses Using real-time PCR (TaqMan®) in Varroa 



PLOS ONE I www.plosone.org 



15 



July 2014 I Volume 9 | Issue 7 | e103073 



Honeybee Colony Disorder in Crop Areas 



Mites Collected from a Thai Honeybee (Apis mellifera) Apiary. Journal of 
Invertebrate Pathology: 91: 69-73. 

46. Martin SJ, Highfield AC, Brettell L, ViUalobos EM, Budge GE, et al. (2012) 
Global honey bee viral landscape altered by a parasitic mite. Science 336: 1304- 
1306. 

47. Zuur AF, leno EN, Elphick CS (2010) A protocol for data exploration to avoid 
common statistical problems. Methods Ecol Evol 1: 3—14. 

48. Beekman M, Ratnieks FLVV (2000) Long-range foraging by the honey-bee. Apis 
mellifera L. Funet Ecol 14: 490-496. 

49. Steffan-dewenter 1, Kiihn A (2003) Honeybee foraging in differentially 
structured landscapes. Proc R Soc Lond B Biol Sci 270: 569-575. 

50. Pinheiro JC. Bates DM (2000) Mixed-Effects Models in S and S-PLUS. 
Springer. 560 p. 

5 1 . Fox J (2002) An R and S-Plus companion to applied regression. Sage Publications, 
Inc. Available: http://books.google.be/booksPhl = en&lr = &id = xWS8kgRjGcAC 
&oi = fnd&pg = PR9&dq = Fox+2002+An+R+and+S+Plus+&ots = o5qCjIMtOb 
&sig = lvi2jQFv4iDIBrvOQtaWiRIlRBM. Accessed 2012 Sep 18. 

52. Bates D (2013) Linear mixed model implementation in lme4. Ms Univ Wis. 
Available: http:/ / dcbianmirror.wwi.dk/cran/web/packages/lme4/vignettes/ 
Implementation.pdf Accessed 2014 Feb 3. 

53. IRM IR dc M (2013) Bilan climatologique annuel. Available: http://ww-\v. 
meteo.be/ mcteo/vicw/fr/ 1 3 1 7239-Bilan-l-climatologiquc+annucl.html. Ac- 
cessed 2013 Dec 15. 

54. De Smet L, Ravoet J, de Miranda JR, Wenseleers T, Mueller MY, et al. (2012) 
BeeDoctor, a versatile MLPA-based diagnostic tool for screening bee viruses. 
PloS One 7: e47953. 

55. Bailey L, Ball BV, Perry JN (1981) The prevalence of viruses of honey bees in 
Britain. Ann Appl Biol 97: 109-118. 

56. Verhoeven B (2014) Projet pilote de surveillance de la mortalite des abcillcs 
2012—2013 Premiers resultats. AFSCA. Journe apicole 26 Janvier 2014, Namur. 

57. Chauzat M-P, Faucon J-P. Martel A-C, Lachaize J, Cougoule N, et al. (2006) A 
survey of pesticide residues in pollen loads collected by honey bees in France. 
J Econ Entomol 99: 253-262. 

58. Chauzat M-P, Faucon J-P (2007) Pesticide residues in beeswax samples collected 
from honey bee colonies (Apis mellifera L.) in France. Pest Manag Sci 63: 1 100- 
1106. doi:10.1002/ps.l451. 

59. Wiest L, Bulete A, Giroud B, Fratta C, Amic S, et al. (2011) Multi-residue 
analysis of 80 environmental contaminants in honeys, honeybees and pollens by 
one extraction procedure followed by liquid and gas chromatography coupled 
with mass spectrometric detection. J Chromatogr A 1218: 5743-5756. 
doi:10.1016/j.chroma.201 1.06.079. 

60. Johnson RM, Pollock HS, Berenbaum MR (2009) Synergistic interactions 
between in-hive miticides in Apis mellifera. J Econ Entomol 102: 474—479. 

61. Yu SJ (2008) The toxicology and biochemistry of insecticides. CRC N Y NY 
USA. Available: http:/ /www.lavoisier.fr/hvre/notice.asp?ouvrage = 1795057. 
Accessed 2014 Jan 31. 

62. DcGrandi-Hoffman G, Chen Y, Simonds R (2013) The Effects of Pesticides on 
Queen Rearing and Virus Titers in Honey Bees (Apis mellifera L.). Insects 4: 
7 1-89. doi: 1 0.3390/insects40 1007 1 . 

63. Pesticide Properties Database (PPDB) (2014) Boscalid (Ref: BA8 51 OF). 
Available: http://sitem.herts.ac.uk/aeru/ppdb/en/86.htm. Accessed 2014Jan 
3. 

64. Aubee C, Lieu D (2010) Environmental Fate and Ecological Risk Assessment for 
Boscahd New Use on Rapeseed, Including Canola (Seed Treatment). U.S. EPA - 
Washington. 

65. Suchail S, Gucz D, Belzunccs LP (2001) Discrepancy between acute and chronic 
toxicity induced by imidacloprid and its metabolites in Apis mellifera. Environ 
Toxicol Chem SETAC 20: 2482-2486. 

66. Belzunccs LP, Colin ME (1992) Les phenomenes d'alfaiblissement de ruchers. 

67. Meled M, Thras^'voulou A, Belzunccs LP (1998) Seasonal variations in 
susceptibility of Apis mellifera to the synergistic action of prochloraz and 
deltamethrin. Environ Toxicol Chem 17: 2517-2520. 



68. Thompson HM (2003) Beha\'ioural elfcets of pesticides in bees— their potential 
for use in risk assessment. Ecotoxicology 12: 317—330. 

69. Mussen EC, Lopez JE, Peng CY (2004) Effects of selected fungicides on growth 
and development of larval honey bees, Apis mellifera L.(Hymenoptcra: Apidae). 
Environ Entomol 33: 1151-1154. 

70. Zhu W, Schmehl DR, MuUin CA, Frazier JL (2014) Four Common Pesticides, 
Their Mixtures and a Formulation Soh'cnt in the Hive Environment Have High 
Oral Toxicity to Honey Bcr Lar\ar. PloS One 9: e77547. 

71. Batra LR, Batra SWT, Bohari CiL ^1973) The mvcoflora oi' domesticated and 
wild bees (Apoidra). AKeopathol Mycol Appl 49: 13-44. 

72. YoderJA, Jajack AJ, Rosselot AE, Smith TJ, Yerkc MC, ct al. (2013) Fungicide 
Contamination Reduces Bcnelicial f ungi in Bee Bread Beised on an Area-Wide 
Field Study in Honey Bee, Apis um/lijem, Colonies. J Toxicol Environ Health A 
76: 587-600. doi:10.1080/15287394.2013. 798846. 

73. Anderson KE, Sheehan IH, Eckhokn BJ, Mott BM, DeGrandi-Holiman G 
(201 1) An emerging paradigm of colony health: microbial balance of the honey 
bee and hive (Apis mellifera). Insectes Sociaux 58: 431—444. doi:10.1007/ 
S00040-0 11-0194-6. 

74. Sartor RB (2008) Therapeutic correction of bacterial dysbiosis discovered by 
molecular techniques. Proc Nad Acad Sci 105: 16413-16414. 

75. Di Pasquale G, Salignon M, Le Conte Y, Belzunccs LP, Decourtye A, et al. 
(2013) Influence of Pollen Nutrition on Honey Bee Health: Do Pollen Quality 
and Diversity Matter? PLoS ONE 8: c72016. doi:10.1371/journal.ponc. 
0072016. 

76. Brodschneider R, Crailsheim K (2010) Nutrition and health in honey bees. 
Apidologie 41: 278-294. 

77. Standifcr LN (1967) A comparison of the protein quaUty of pollens for growth- 
stimulation of the hypopharyngeal glands and longevity of honey bees, Apis 
mellifera L. (Hymenoptera: Apidae). Insectes Sociaux 14: 415-425. 
doi:10.1007/BF02223687. 

78. Liming W, Jinhui Z, Xiaofeng X, Yi L, Jing Z (2009) Fast determination of 26 
amino acids and their content changes in royal jelly during storage using ultra- 
performance liquid chromatography. J Food Compos Anal 22: 242-249. 

79. Kubik M, Nowacki J, Pidek A, Warakomska Z, Miehalczuk L, et al. (1999) 
Pesticide residues in bee products collected from cherry trees protected during 
blooming period with contact and systemic fungicides. Apidologie 30: 12. 
doi: 1 0. 1 05 1 /apido: 1 9990607 . 

80. Kubik M, Nowacki J, Pidek A, Warakomska Z, Miehalczuk L, et al. (2000) 
Residues of captan (contact) and difenoconazole (systemic) fungicides in bee 
products from an apple orchard. Apidologie 31: 531-542. 

81. Pesticide Properties Database (PPDB) (2014) Pyrimethanil (Ref: SN 100309). 
Available: http://sitem.herts.ac.uk/aeru/ppdb/en/573.htm. Accessed 2014Jan 
3. 

82. Ravoet J, Maharramov J, Mecus 1, De Smet L, Wenseleers T, ct al. (2013) 
Comprehensive bee pathogen screening in Belgium reveals Crithidia mellificac 
as a new contributory factor to winter mortality. PloS One 8: e72443. 

83. Oomen PA, Belzunccs L, Pelissier C, Lewis G (2001) Honey bee poisoning 
incidents over the last ten years, as reported by bee keepers in The Netherlands. 
Les Colloques 98: 129-136. 

84. Brittain CA, Vighi M, Bommarco R, SerteleJ, Potts SG (2010) Impacts of a 
pesticide on pollinator species richness at different spatial scales. Basic Appl Ecol 
11: 106-115. 

85. CRA-API (201 1) Efiects of coated maize seed on honcv bees. Report based on 
results obtained from the first year of activity of the APENET project. Italy. 

86. Maxim L, van der Sluijs J (2013) Seed-dressing systemic insecticides and 
honeybees. Late Lessons from Early Warnings: Science, Precaution, Innovation. 
Copenhagen. 401-438. Available: http:/ /www.beekeeping.com/articles/us/ 
late_lessons_from_early_warnings_2.pdf. Accessed 2014 Feb 10. 

87. Holzschuh A, Steffan-Dewentcr I, Tscharntke T (2008) Agricultural landscapes 
with organic crops support higher pollinator diversity. Oikos 117: 354—361. 

88. Kleijn D, Kohler F, Baldi A, Batary P, Concepcion ED, et al. (2009) On die 
relationship between farmland biodiversity and land-use intensity in Europe. 
Proc R Soc B Biol Sci 276: 903-909. 



PLOS ONE I www.plosone.org 



16 



July 2014 I Volume 9 | Issue 7 | e103073 



